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^ ■ ABSTRACT 
C/3 . 

. We present medium resolution spectropolarimetry and long term photo- 

polarimetry of two massive post-red supergiants, IRC +10420 and HD 179821. The 
data provide new information on their circumstellar material as well as their evolution. 
^ ' In IRC +10420, the polarization of the Ha line is different to that of the continuum, 

00 ' which indicates that the electron-scattering region is not spherically symmetric. The 

, observed long term changes in the polarimetry can be associated with an axi-symmetric 

OO ■ structure, along the short axis of the extended reflection nebulosity. Long term pho- 

I tometry reveals that the star increased in temperature until the mid-nineties, after 

which the photospheric flux in the optical levelled off. As the photometric changes are 
' mostly probed in the red, they do not trace high stellar temperatures sensitively. And 

OO ! so, it is not obvious whether the star has halted its increase in temperature or not. For 

■ HD 179821 we find no polarization effects across any absorption or emission lines, but 

observe very large polarization changes of order 5% over 15 years. During the same 
period, the optical photometry displayed modest variability at the 0.2 magnitude level. 
' This is unexpected, because large polarization changes are generally accompanied by 

I strong photometric changes. Several explanations for this puzzling fact are discussed. 

Most of which, involving asymmetries in the circumstellar material, seem to fail as 
there is no evidence for the presence of hot, dusty material close to the star. A caveat 
is that the sparsely available near-infrared photometry could have missed periods of 
strong polarization activity. Alternatively, the variations can be explained by the pres- 
ence of a non-radially pulsating photosphere. Changes in the photometry hint at an 
increase in temperature corresponding to a change through two spectral subclasses 
over the past ten years. 

Key words: Techniques: polarimetric - Stars: evolution - Stars: circumstellar envi- 
ronment - Stars individual : IRC +10420, HD 179821 



1 INTRODUCTION 

Due to the steepness of the Initial Mass Function, massive 
stars (> 8 M©) are extremely rare. This is exacerbated by 
their comparatively short lifetimes. Yet, although rare, these 
objects have a crucial impact on the interstellar medium 
due to their strong winds and high mass-loss rates, and can 
dominate the light output of entire galaxies. 

An area of current interest is that massive evolved 



stars are often surrounded by bi-polar nebulae ('Weid l20"03l ). 
Later in the evolution of a star, using spectropolarime- 
try, it is now established that the eje cta of supe r novae 
deviate from spher ical symmetry (e.g. IWang et al. I l2003l : 
[Leonard et all |2005| ). It is as yet unclear whether this is 
due to asymmetric explosions, axi-symmetric stellar winds 
or a pre-e xisting density contrast in the surrounding ma- 
terial (e.g. iDwarkadas fc Balic3ll998l : lDwarkadas fc Owockil 
l2002f) . Wind-axisymmetry may imply fast rotation, and be 
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related to the beaming of SN explosions, which may be the 
origi n of the extremely luminous, beamed gamma-ray bursts 
(e.g. iMeszarodlioO^ : iMazzah et al]l2003l ). 

Here, we address the issue by investigating the circum- 
stellar ejecta of two yellow hypergiants, IRC +10420 and 
HD 179821. These objects are thought to have evolved off 
the post-Red Supergiant branch and are still surrounded by 
mass ejected during a previous mass losing phase. O nly a few 
yellow hypergiants are known (see the review by Ide Jaeen 
[l998), and the number of such hypergiants with circumstel- 
lar dust is even smaller - only IRC +10420 and HD 17982 1 
belong to this class (see for example lOudmaiier et allbOOSl ). 
Therefore, the study of these two unique objects is impor- 
tant in its own right. 

As these stars are distant (3-5 kpc), the direct imag- 
ing of their innermost regions is currently beyond the 
reaches of current technology, although interfero metry is 
start ing to resolve the winds of evolved stars (do Wit et al.l 
I2OO8I ). Observing with spectropolarimetry allows us to probe 
regions much closer to the star still. Spectropolarime- 
try was first effectively used in the study of classical Be 
stars using the presence of 'l ine effects' (iPoeckerd 1 19751 : 
IPoeckert fc Marlboroughl Il976l ). These are changes in po- 
larization across spectral lines that have an emission com- 
ponent. They occur because emission-line photons arise over 
a larger volume than the stellar continuum photons. Conse- 
quently, the emission-line photons undergo fewer scatterings 
as they 'see' fewer electrons, resulting in a lower polariza- 
tion than the continuum. We normally only observe a net 
polarization change if the geometry of the electron scatter- 
ing region is aspherical. Many authors have confirmed that 
this technique provides evidence that envelope geometries 
around Be stars are indeed disk-like (|Doughertv fc Tavlorl 
1 19921 : iQuirrenbach et al.lll997l ; IWood et al.lll997l ). More re- 
cently the technique has been used to investigate the geom- 
etry of cir cumstellar material around H erbig Ae/Be stars. 
Studies bv lOudmaiier fc Drewl (|l999l ) and lVink et all (|2002l ) 
show that most Herbig stars exhibit line effects, indicating 
asp herical elec t ron-sc attering regions. For a recent review, 
see lOudmaiierl (|2007l ). With regard to evolved stars, Davies 
et al. (2005) conducted a study of Luminous Blue Variables 
(LBVs) using spectropolarimetry. They found that 50% of 
the objects observed exhibited polarization changes across 
Ha, indicating that some asphericity lies at the base of the 
stellar wind. Furthermore, they found several objects for 
which the position angle varied randomly with time, leading 
them to conclude that the wind around these stars is clumpy 
(see also iNordsieck et al.|[200lh . 

IRC +1 0420 is now well accepted as a massive, evo lved 



object (e.g. Jones et al.l 1 19931 : lOudmaiier et aL 1 1 19961 and 



iHumphrevs et aL I l2002t ). This is mainly based on its large 
distance, high outflow velocity (40 kms~^) and high lumi- 
nosity implied from the hypergiant spectrum. The situation 
for HD 179821 is less certain. The presence of non-radial 
pulsations coupled with comparati vely modest photometri c 
changes suggest a massive nature l|Le Coroller et aLll2003l ). 
Furthermore, its circumstellar material has a large expan- 
sion velocity of 30 km s~^, as measured in CO, su ggesting 
the star is a supergiant (|Kastner fc Weintraubl 19951 ). On the 
other hand, the overabundance of s-process elements and the 
low metallicity suggest HD 179821 is perhaps a lower mass 



post-AGB star JZacs et al.l 1 19961 : iReddv fc HrivnakI Il999l : 
iThevenin et aLlbOOd ). 

Although the latter's nature is a bit more uncertain, 
there are some striking similarities between IRC +10420 
and HD 179821. When observed as part of a larger sample of 
post-AGB stars, these two objects are often markedly differ- 
ent from the rest. In particular, their high outflow velocities 
(the average outflow velocity for post-AGB and AGB stars 
is 15 kms~^) require much higher luminositie s if powered 
by radiation pressure alone (|Habing et al.|[l993 ). They were 
the only objects t hat showed extensive reflectio n nebulae in 
a larg e survey by iKastner fc Weintraubl l|l995l ) . lJura et al] 
l|200ll ) point out the enormous difference between the space 
velocities of IRC +10420 and HD 179821 when compared 
against low mass post-AGB stars. Furthermore, both objects 
have an exceptionally strong Ol 777 4 tr iplet absorptio n fea- 
ture indicati ng a high lumin osity (|Hurnphreys et al] 1973 
and [keddv fc HrivnakI Il999l . based on ISlowik fc PetersonI 
Il995h . We therefore proceed with both objects and assume 
they are evolved post-Red Supergiants. 

Recently, both objects have been observed at arcsecond 
resolution in CO bv .Castro-Carrizo et a l. (2007) who found, 
in accordance with previous estimates, that their mass loss 
rates exceeded 10~* MQyr~^ when they were in the RSG 
phase. The envelopes show mild deviations from spherical 
symmetry in their data. 

This paper is organised as follows. In Section 2, we re- 
view the experimental setup and explain how the data has 
been reduced. We present our results for each object in Sec- 
tion 3, and use the new spectropolarimetric data together 
with past polarization measurements to investigate the na- 
ture of the circumstellar environments around each of the 
stars, which is discussed in Section 4. We conclude in Sec- 
tion 5. 



2 OBSERVATIONS 

We describe our most recent data in detail. Similar observa- 
tions taken earlier are briefly discussed, with relevant differ- 
ences highlighted. The linear spectropolarimetric data were 
taken on the night of 30 September 2004 using the ISIS spec- 
trograph on the 4.2m William Herschel Telescope (WHT), 
La Palma. A MARCONI2 CCD detector with a R1200R 
grating was used. This yielded a spectral coverage of 6150- 
6815 A and gave a spectral resolution of 34 kms~^ at Hq. 
The seeing was about 2" , and a slit with a width of 1 arcsec 
was used. The linear polarimetric component of the data was 
analysed using the polarization optics equipment present on 
the ISIS spectrograph. The object and the sky background 
were simultaneously observed using additional holes in the 
dekker mask. A calcite block was then used to split the rays 
into two perpendicularly polarized beams (the o and e rays). 
A complete data set therefore, consists of four spectra ob- 
served at 0° and 45° (to measure Stokes Q) and 22.5° and 
67.5° (to measure Stokes U). 

All data reduction steps used the FIGARO software main- 
tained by Starlink and included bias subtraction, cosmic ray 
removal, bad pixel correction, spectrum straightening and 
flat fielding. Wavelength calibration was carried out using 
observations of a Copper- Argon lamp taken throughout the 
run. The data were then imported into the package CCd2pol 
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Object Telescope Date Julian Date %P P.A. (Deg) %P (Ha) P.A. (Ha) (Deg) 

HD 179821 



IRC +10420 



AAT 


15/09/02 


2452533 


1.99 ± 0.11 


40 ± 2 






WHT 


30/09/04 


2453279 


2.00 ± 0.15 


36 ± 3 






NOT 


16/01/98 


2450830 


1.95 it 0.05 


174 d= 1 






NOT 


16/05/98 


2450950 


1.80 ± 0.03 


174 ± 1 






AAT 


18/09/02 


2452536 


2.12 ± 0.11 


173 ± 2 


1.28 ± 0.11 


10 ± 2 


AAT 


15/08/03 


2452867 


2.35 ± 0.11 


179 ± 1 


1.61 ± 0.11 


11 ± 1 


WHT 


30/09/04 


2453279 


3.40 ± 0.15 


174 ± 2 


1.93 ± 0.15 


11 ± 2 



Table 1. New polarimetric observations of HD 179821 and IRC +10420 measured in the R band. For the spectropolarimetric data 

(those taken at the AAT and WHT), the polarization was measured in the continuum region close to Ha, while the final columns are 
the polarizations at the line centers. The Ha line-center polarization of HD 179821 was not calculated as no line effect was observed, and 
therefore this measurement would be equal to the continuum polarization. The systematic error of the AAT and WHT spectropolarimetric 
data is estimated to be of order 0.10% and 0.15%, respectively. 



JulicLii Date 


Telescope 


rS — y 


1/ 
V 


V — rt 


I? T 


T 

J 


ZJ 
±1 


I\ 


2450268 


CST 










5.62 






2450303 


CST 










5.36 


4.40 


3.45 


2450643 


CST 










5.40 


AAA 


3.50 


2450690 


CST 










5.40 


4.45 


3.49 


2450707 


CST 










5.37 


4.43 


3.38 


2450950 


NOT 


2.76 


11.06 


2.40 


1.70 








2450985 


CST 










5.38 


AAA 


3.48 


2451037 


TSAO 


2.58 


11.12 


2.42 


1.60 








2451039 


TSAO 


2.67 


11.11 


2.41 


1.58 








2451042 


TSAO 


2.58 


11.13 


2.42 


1.54 








2451043 


TSAO 


2.73 


11.06 


2.49 


1.61 








2451047 


TSAO 


2.51 


11.01 


2.41 


1.58 


5.63 


4.55 


3.73 


2451048 


TSAO 


2.75 


11.03 


2.47 


1.60 


5.43 


4.57 


3.58 


2451050 


TSAO 


2.67 


11.01 


2.46 


1.66 








2451052 


TSAO 


2.75 


11.03 


2.47 


1.59 


5.35 


AAA 


3.53 


2451057 


TSAO 


2.76 


10.96 


2.47 


1.64 








2451063 


TSAO 


2.60 


10.95 


2.43 


1.59 








2451075 


TSAO 


2.70 


11.16 


2.52 


1.60 








2451082 


TSAO 


2.66 


11.09 


2.47 


1.59 








2451083 


TSAO 


2.72 


11.05 


2.48 


1.60 








2451087 


CST 










5.24 


4.29 


3.33 


2451099 


TSAO 


2.62 


11.16 


2.56 


1.66 








2451100 


TSAO 


2.71 


11.16 


2.54 


1.66 








2451103 


TSAO 


2.76 


11.13 


2.51 


1.54 








2451104 


TSAO 


2.72 


11.03 


2.48 


1.59 








2451147 


CST 










5.35 


4.42 


3.46 


2451245 


TSAO 










5.32 


4.47 


3.56 


2451292 


CST 










5.40 


4.45 


3.57 



Table 2. Photometry of IRC +10420. The data come from the Carlos Sanchez Telescope (CST), the Nordic Optical Telescope (NOT) 
and the Tien-Shan Astronomical Observatory (TSAO), see text for details. Typically, the photometric errors are of order 0.01-0.03 
magnitude. 

a detection of 10^ photons, and significantly smaller when 
large parts of the spectra arc binned. However, the presence 
of systematic errors such as scattered light and instrumental 
polarization often outweigh these photon-noise errors and we 
estimate the errors in the present data to be of order 0.15%. 
The observations of two zero-polarization objects show that 
the instrumental polarization is less than 0.15%. We have 
made no attempt to correct for either instrumental or inter- 
stellar polarization (ISP), as these only add a constant (Q,U) 
vector to the data and therefore will not affect whether we 
observe a line effect. 



(also maintained by Starlink), to produce the Stokes Q and 
U parameters. The degree of polarization and its PA can 

then be found from 

= 0" + (1) 

e = 0.5 X arctan(^) (2) 

Finally, the data was position angle-calibrated using a 
set of polarization standard stars observed during the run. 

Polarimetric data should only be limited by photon 
statistics. For example we expect an error of about 0.1% for 
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IRC + 10420 09/02 error = 0.1 IRC+10420 08/03 error = 0,1 IRC+10420 09/04 error = 0.3 




6520 6540 6560 65B0 6600 6520 6540 6560 6580 6600 6520 6540 6560 6580 6600 

Wavelength (A) Wavelength (A) 'A'avelength (A) 




1 1.5 2 1.5 2 2. S 12 3 4 

q(%) q(%) q(%) 

Figure 1. Spcctropolarimctry of IRC +10420 around Ha at different epoclis. Tlie upper diagrams are 'triplets' of the spectropolarimctric 
data, which were taken one year apart from one another. The bottom panel of the triplot presents the direct intensity spectrum, the data 
normally obtained from spectroscopy. The middle panel shows the polarization measured as a function of wavelength. The upper panel 
shows the corresponding polarization angle. The data are adaptively binned such that each bin has the same observational error derived 
from Poisson statistics. Here we have chosen errors of 0.1%, 0.1% and 0.3% respectively for the 2002, 2003, and 2004 data. Below the 
triplets, the Stokes Q,U vectors are plotted against each other, and binned to the same precision as the spectra. Depolarization around 
Ha is seen each time and the associated excursions in Q-U space all point in the same direction. 



2.1 Additional Data 

2.1.1 Spectropolarimetry 

Further data of IRC +10420 and HD 179821 was obtained 
during previous observing runs. Spectropolarimetric obser- 
vations of IRC +10420 were obtained using the 3.9m Anglo- 
Australian Telescope (AAT) on 18/9/2002 and 15/8/2003. 
HD 179821 was observed on 15/9/2002 at the AAT. The 
2002 data are taken with a lower spectral resolution (2.1 A) 
and the 2003 data are taken at comparable resolution to the 
current WHT data. We refer to FDavies et al. (2005. ) for more 
details on the observations taken at the AAT in 2002 and 
2003. 



servations are described in lOudmaiier et all (|200ll ). John- 
son BVRIJHK observations were obtained in 1998 on the 
1 m telescope at the Tien-Shan Astronomical Observatory 
(Kazakhstan) with the 2d-channel p hotometer-pola r imete r 
FP3U of the Pulkovo Observatory (iBergner et al.1 119881 ). 
Further near-infrared photometry was obtained employing 
the 1.5m Carlos Sanchez Telescope in Tenerife. Details of 
the ob serving procedures can be found in lKerschbaum et al.l 
(l2006t ). The results are presented in Tables [T] and H) R band 
polarimetric data and optical and near-IR photometric data 
taken from the literature, for both objects, is described in 
Tables [3] and 3] respectively. 



2.1.2 Photometry and polarimetry 

We also obtained additional photometric and polarimetric 
data of IRC +10420. Hitherto unpublished broad band po- 
larimetry was obtained in January and May 1998 using the 
2.5m Nordic Optical Telescope in La Palma, and on one 
of these nights B VRI photometry were obtained. These ob- 



3 RESULTS 
3.1 IRC +10420 

3.1.1 Spectropolarimetry 

In Fig. [l]we have plotted the polarimetric data around Ha 
taken in 2002, 2003 and 2004. These data show 'triplots' 
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Object 


Date 


N 


Source 


Method 


HD 179821 












May 1989 - May 1991 


4 


Pajrthasajrathy et al. (2005) 


Broadband 




October 1991 


1 


Trammell et al. (1994) 


Spectropolarimetry 




August 1993 


1 


HPOL database, Johnson (priv. communication) 


Spectropolarimetry 




June 1997 - October 1998 


30 


Melikian et al. (2000) 


Broadband 


IRC +10420 












May 1976 - June 1976 


1 


Craine ct al. (1976) 


Broadband 




May 1989 


1 


Johnson & Jones (1991) 


Broadband 




October 1991 


1 


Jones ct al. (199.3) 


Broadband 




October 1991 


1 


Trammell et al. (1994) 


Spectropolarimetry 



Table 3. R, band polarimetric observations of HD 179821 and IRC +10420 taken from the literature. The observational period is given 
in the second column. The last three columns show the number of observations taken, the literature source of the data and the method. 
Although Craine et al. (1976) made several observations of IRC +10420, these values were averaged to a single data point in our study. 
We measured the continuum close to Ha from the spectropolarimetric data in order to facilitate comparison with the R band data. 



Object 


Date 


Band 


N 


Source 


HD 179821 












March 1988 - April 1988 


V,J,K 


1,2 


Hrivnak et al. (1993) 




May 1990 - October 1999 


V 


169 


Arkhipova et al. (2001) 




May 1990 - November 2005 


V 


254 


ASAS catalogue 




June 1992 


J,K 


1 


Kastner & Weintraub (1995) 




October 1998 


J,K 


1 


2MASS database 


IRC +10420 












May 1996 - June 1996 


J,K 


2,2 


Humphreys et al. (1997) 




April 2001 


J,K 


4,5 


Kimeswenger et al. (2004) 




April 2002 - October 2004 


V 


104 


ASAS Catalogue 



Table 4. A summary of recent photometric observations of IRC +10420 and HD 179821. In the fourth column we give the number of 
observations taken in each passband, respectively. For IRC +10420, this acts as an update to the previous tables compiled by Jones et 
al. (1993) and Oudmaijer et al. (1996). 
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Figure 2. Left : The polarization and position angle vajriability of IRC +10420 in the iJ-band. Both the polarization and PA show 
an increase with time. The bottom 2 panels show the evolution of the Stokes Q and U vectors. Right : An updated version from the 
photometric variations presented in Oudmaijer et al (1996). The change in the 7 band, which traces the stellar photosphere, indicates a 
temperature increase of the star. Also included is the evolution of the Ha Equivalent Width which displays a marked strengthening over 
the years. 
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Figure 3. A Q-U diagram of IRC +10420. Tlie fi-band (Q,U) 
data-points come from the literature and our new broadband po- 
larimetry. The dashed line represents the vec tor towards the ISP , 
denoted by the crossed circle, as derived by Ijones et al 
The continuum values from the spectropolarimetric data are rep- 
resented by dots and the straight lines show the excursions to the 
line center polarization measured from the spectropolarimetry. 

of the spectropolarimetry, along with graphs plotting the 
Stokes Q,U vectors against each other. 

IRC -1-10420 displays Ha emission with a total equiva- 
lent width (EW) of -50 A. The emission line is unresolved in 
the lower resolution 2002 data. The other, higher resolution, 
data indicate a double-peaked and variable emission-line 
profile. In 2003, the red peak is stronger than the blue peak, 
but in 2004 the blue peak has become the stronger, while 
the central absorption component of the double-peaked line 
has become deeper. Comparison with published data reveals 
that the red peak is very strong in the present spectrum. 
lOudmaiierl (ll998D show data taken in 1993 of the object, 
where the blue peak is 2 to 3 times stronger than the red 
one. The EW of the line was -25A, about half that observed 
now. This evolving blu e-to-red peak ratio was also noted by 
iKlochkova et all |2002| ). 

The spectropolarimetric data show large polarization 
changes across the line. In each case, this line effect consists 
of a depolarization over the central part of the line coupled 
with a rotation to larger angles. The associated Q-U plots 
(Fig. [T} also clearly show the line effects. The cluster of data 
points in the lower right hand corners of the graphs are due 
to the stellar continuum, while the excursions to the upper 
left hand corners trace the polarization over the Ha line. 
At first glance, the excursions might seem to suggest the 
presence of structure. However, given that the binning is at 
the 0.1% and 0.3% level respectively, and that the 'structure' 
is mostly due to single points offset by less than 3(j from 
the global trends, we proceed under the assumption that 
the excursions are linear. The intrinsic polarization angle 
for IRC -f 10420 is measured by the excursion from the line- 
center (which has a lower contribution by the polarization 
due to electron scattering) to the continuum. Using Eq. 2 
this is 158°, with an estimated error of 2°, and it does not 
seem to change with time. The magnitude of the excursion 



gives an indication of the strength of the polarization change 
and is similar in 2002 and 2003, but almost twice as large in 
2004. Strictly speaking the length of the vector represents a 
lower limit to the real depolarization because a finite spectral 
resolution can wash out the line-effect. However, the 2002, 
lower resolution data has a depolarization vector of similar 
length to the higher-resolution, 2003 data. It is therefore 
plausible that the depolarization (and thus the polarization 
of the continuum due to electron scattering) is significantly 
stronger in the 2004 data. 

The continuum polarization towards IRC -1-10420 con- 
sists of three components. Firstly, there is interstellar polar- 
ization due to dichroic absorption of interstellar dust. This 
is constant with time and is constant over a small part of 
the spectrum. In data such as in Fig. [T] this would add a 
constant vector to the intrinsic Q,U spectrum. Importantly, 
the shape of the QU behaviour is not affected by the ISP. 
Secondly, some polarization may be due to scattered hght 
reflected off circumstellar dust grains. The presence of scat- 
tered light is revealed by iKastner fc Weintraubl l|l995h who 
show that the polarization is significant at large distances 
from the star. Later HST images show substantial extended 
emission in blue light which is interpr eted as s cattered light 
as we ll l|Humphrevs et aLlll997l . |2002| : see also iDavies et al.l 
l2007f ). When obtaining spatially unresolved data, such as 
here, any net polarization is due to the asymmetry of the 
scattering material on the sky. Thirdly, the continuum light 
can be scattered off free electrons that are very close to the 
star. Here, imaging polarimetry can not achieve the spatial 
resolution required to study the electron scattering region. 
The line-effect however readily reveals its presence. 

3.1.2 Long term photo-polarimetry 

Our new data are combined with literature data and are 
shown in Figure [5] Over the past 30 years, both the po- 
larization and the angle have increased. To illustrate this, 
the Q and U Stokes vectors are also plotted. The Stokes Q 
parameter increases by about 2%, while the U vector essen- 
tially remains constan t (see also Fig. Here we note the 
iTrammell et all l| 19941 ) data. Their data point (3 . 7% at 149°) 
was taken in the same month as the I Jones et al.l l| 19931 ) data 
which is 2.0% at 174°. The Trammell et al. measurement 
strongly deviates from the general, 30-year, trend in P, O 
and U and is also significantly different from data taken in 
the same month. We will proceed assuming that the trend 
of increasing polarization with time reflects the true long 
term polarization behaviour of the object, and have a caveat 
that, perhaps, the evolution is more complex. In Fig. |3]we 
plot the broad-band Stoke s Q and U vectors against one 
another. iJones et al.l (|l993l ) derive an interstellar R band 
polarization of 8 per cent towards IRC -1-10420 - a value 
that is not inconsistent with the large reddening towards 
the object. From their figures we measure their ISP to be 
Qr,isp ~ 7.5%, Urjsp ^ 2.7%. This value is indicated in 
Fig. |3l Assuming the Ha line center is not polarized due to 
electron scattering, its net polarization is due to a combina- 
tion of ISP and scattering off the circumstellar dust. We can 
thus get a (very) rough estimate of the circumstellar dust 
polarization by taking the difference between the Ha line 
center polarization and the ISP value. Our observed line 
center polarizations are close each other. Remarkably, the 
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recent Hq line polarization values all lie on the vector con- 
necting the ISP with the origin. Their average distance to 
the ISP value is {Q = -6.0%, U = -2.15%), corresponding 
to a circumstellar dust polarization of {P = 6.4%, — 9.9°). 
The true value may have a larger polarization at a smaller 
angle, as the polarization in the center of the line is proba- 
bly affected by the resolution. Note that an interesting sit- 
uation has arisen where the circumstellar dust polarization 
effectively cancels part of the ISP resulting in smaller net 
polarization values. Using Eq. 2, we measure an intrinsic po- 
larization angle for the variations of 6 ±2°. This is very close 
to the orientation found for the circumstellar dust polariza- 
tion described above and pinpoints the changes in circum- 
stellar dust polarization as a likely cause for the observed 
variations. 

The photometric variability has bee n discussed exten- 
sively bv ljones et al.l (1993) and Oudm aiier et al.l l|l996l ). 
Here we give an update on the long term variability (see 
Table 3}. The K band data show a steady fading starting 
in the mid-seventies. This dimming is mimicked by the data 
in the J band, w hich is known to trace the stellar photo- 
sphere lOudinaii er et al.lll996l ). After around 1995, the J 
band magnitudes appear to reach a steady state. 

Long term changes have also b een observed in the hy - 
drogen recombination lines. In 1992, iQudmaiier et al.l (|l994l ') 
discovered NIR recombination emis sion lines that were pre - 
viously seen in absorption in 1984. Ilrvine fc Herbi3 (|l986l ) 
detected Ha emission in 1986, while previous observations 
did not show emission. Therefore, the onset of hydrogen re- 
combination e mission can be traced to a period between 
1984 and 1986 l|Oudmaiier et al.' 'l994Y This coincides pre- 
cisely with a period of strong photometric changes. 

EW data of Ha were collected from the literature as 
well. It is not trivial to obtain these data because of the 
large width of the line. For example, EW measurements 
from high resolution echelle data are not possible as the 
free spectral range is too small to properly define a base- 
line and measure the equivalent width accurately. We there- 
fore restrict ourselves to results from spectra with a large 
free spectral ra nge - essentially no echel le observations. Dat a 
are taken from |Oudmaiier et al.1 (|1994| ) ■ lOudmaiieij ll 19981) . 
iHumphrevs et aL ( 20021 ). this paper and lDavies et al.l ( 20071 ) 
respectively. The results are also shown in Fig. [51 Ha has 
been increasing in strength, going from upper limits of the 
EW in the seventies to -23 A in the early nineties to -50 A 
in 2005. As the optical brightness has remained constant in 
that time, these EW changes refiect true line flux changes 
and imply that more ionized material is present. 



3.2 HD 179821 

3.2.1 Spectropolarimetry 

Spectropolarimetric data of HD 179821 are shown in Fig. 3] 
No polarization changes are visible across the faint Ha line. 
The 2002 data has a higher signal-to-noise, but a weaker 
Ha emission than in 2004. There are few reports of the Ha 
line in the literature, but it would appear it is variable in 
the sense that it varies between weak emission above the 
continuum (i.e. 10-20% emi ssion above the continuum was 
observed bv lZacs et al.|[l996l ) and weak emission that hardly 
fills in the underlying absorption as observed here. 



The weak emission implies that the number of free scat- 
tering electrons is rather limited, which is not inconsistent 
with the low surface temperature of the star (its spectral 
type is G4 0-Ia. ,Keenan fc Mc Neil 1989). Therefore any con- 
tinuum polarization due to electron scattering is bound to 
be low. In such a case, extremely high SNR data are re- 
quired to reveal a line depolarization, if it were present at 
all. Alternatively, the hydrogen recombination line emission 
in HD 179821 could be due to photospheric shocks. The ob- 
ject has been found to s how periodic photomet ric variations 
of 140 and 200 days (see lLe Coroller et al.|[2003h . These vari- 
ations are most likely due to pulsations, and it is common 
to observe Ha emission due to sho cked layers in the ph oto- 
sphere of cool, pulsating stars (e.g. ISchmidt et al.ll2004l and 
references therein). 

3.2.2 Long term photo-polarimetry 

The photometry and polarimetry from this paper and the lit- 
erature of HD 179821 are plotted i n Fig. [5l The longer term 
optical variability is discussed by lArkhipova et al.l (|200lh . 
From 1989 to about 1999 the star became fainter by about 
0.1 magnitude, and t h en re turned to its original magni- 
tude. lArkhipova et al.1 (I2OOII ) suggest the photometry may 
be cyclic on long timescales. They also find that the object is 
bluer in U — B when the star is brightest in V, which they 
attribute to pulsations. This is confirmed by the detailed 
study on shorter term va riations in the V band presented by 
iLe Coroller et al.1 (|2003| ). which are due to pulsations with 
periods of order hundreds of days. Here we report the hith- 
erto unnoticed large change in the J and K bands. These 
became fainter by around 0.4 magnitude in the period 1988 
- 2000. 

Figure [S] shows that HD179821 experiences polarization 
changes (the star is even observed to exhibit no polarization 
at all for four consecutive days in November 1997) , suggest- 
ing that the star is intrinsically polarized. The polarization 
varies more or less randomly between and 2.5 per cent over 
the past 15 years, but we can identify two clear values for 
the PA, a low value around 40° and a high value of around 
120°. When plotting these data in QU space (Fig. (6)) we can 
see that the data points are distributed around what we can 
loosely describe as a straight line through the origin. 

The move through the origin of the QU diagram is re- 
sponsible for the large change in observed polarization an- 
gle. This must have happened at least five times: the first 
occurred between the 1989-1991 IParthasarathv et al.l (|2005l ) 
dataset (PA=40°) and the 1991 Trammcll ct al. (199^ data 
point (PA=120°) taken half a year later, when the object 
moved from the upper half to the lower half in the QU di- 
agram. After this, a move occurred between 1991 and the 
HPOL data taken i n 1993 (John son private communication, 
using HPOL, Wol ff et al.1 1 19961 ) . The object re t urned back 
to the lower half as observed bv lMelikian et al] (|2000l ). who 
also directly detected such an excursion in 1997-1998 (the 
non-detections which lie, by definition, at the origin). Finally 
a crossing happened after their observing run, but before our 
present data was taken. The move through the origin does 
not necessarily refiect a real 90° rotation, as the ISP may 
contribute significantly to the total Q,U vectors observed 
towards the object. 

It is hard to determine the ISP towards any object. 
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Figure 4. As in Fig. [T] but now for HD 179821 in 2002 (left) and 2004 (right). The data has been rebinned such that the Icr error in 
polarization corresponds to the value stated at the top of the triplot, i.e. 0.2% and 0.5% respectively, as calculated from photon statistics. 



but we can make an educated g uess. Since Av ~ 2 for 
HD 179821 ((Hrivnak et all 1 1989 ) and one typically finds 
that the percentaRe polarization is equal t o the extinction 
in magnitudes (e.g. lOudmaiier et al.ll2()()lh . we can expect 
an ISP of up to 2 per cent towards the object. A value for 
the interstellar polarization angle may be estimated from 
the polarization of surrounding field stars. To this end we 
invest igated the compilation of polarized stars due to lHeilej 
The local ISP appears very ordered towards this part 
of the sky. We selected all objects with a polarization mea- 
surement larger than 0% and within a radius of 5 degrees 
from the position of HD 179821. These 22 objects have an 
average polarization angle of 50° with a scatter of 24° . The 
catalogue also lists the photometric distances, and the aver- 
age angle is the same for the 11 furthest (>500pc) objects. 
Whatever the magnitude of the polarization, these values 
indicate that the ISP can be located in the top quadrants of 
Fig.O When taking these angles at face value, and assuming 
the ISP is oriented in the same direction at large distances 
as locally, then the ISP towards HD 179821 is 50° at 2%. 
This is close to the present-day value of HD 179821, and, if 
true, would suggest that HD 179821 is presently in a rela- 
tively quiescent state, yet the polarization is still variable by 
up to a few per cent. The instances when the object occu- 
pies the bottom of the Q U diagram can then be identified as 
periods where the polarization of HD 179821 is particularly 
strong. The optical photometry does not signal any such ac- 
tivity however, and continues its gradual decline to fainter 
magnitudes. 

In short, during the last 15 years, the observed polar- 



ization of HD 179821 crossed the origin in the QU-diagram 
several times. The changes occur along a straight line and 
suggest some preferred orientation, for which we derive a 
very rough estimate for the intrinsic angle of 38°. In addi- 
tion, the polarization changes by 4-5% , while only modest 
changes in the optical photometry of less than 0.2 magnitude 
occur. Near-infrared photometry is sparse, and indicates a 
gradual faintening of the star. The most recent data point 
obtained o n 12 October 1998, wh ich was a few weeks before 
the end of lMehkian et all (|2000l ')'s campaign, i.e. when HD 
179821 was, if the ISP value is correct, still in a high state 
of polarization. 



4 DISCUSSION 

We have observed the only two post-Red Supergiants with 
large infrared excesses using multi-epoch spectropolarime- 
try and supplemented these data with long term photo- 
polarimetry. We find some striking similarities between the 
two objects, but also large differences. We will start with 
the latter. 



4.1 The ionized gas 

The spectropolarimetry reveals a strong line-effect in the 
data of IRC -hl0420, but none for HD 179821. The lat- 
ter displays weak Ha emission and is much cooler than 
IRC -1-10420. It is therefore not surrounded by much ionized 
material. As electron-scattering is the prime mechanism to 
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Figure 5. The variability of HD 179821 with time. Plotted is the 
V-band photometry provided by Hrivnak et al. (1989) [squares], 
Arkhipova et al (2001) [triangles] and the ASAS Catalogue (Poj- 
manski, 2002) [circles]. The error associated with the photometry 
measurements provided by Arkhipova is 0.005 mag and those from 
the ASAS Catalogue are less than 0.06 mag. Plotted in the upper 
two plots are PA and polarization data taken from Parthasarathy 
et al. (2005) [crosses], HPOL database [triangle], Trammell et al. 
(1994) [star], Melikian et al. (2000) [circles] and data presented 
in this paper from 2002 [squares] and 2004 [rhombi]. The PA ex- 
hibits large rotations over a 15 year period, while the polarization 
indicates the presence of short-term variability. The lower panels 
show that HD 179821 has become fainter in J and K. 



produce line-effects, it can then be easily understood why 
there is no line-effect towards HD 179821. 

As described earlier, the intrinsic polarization angle 
due to electron scattering for IRC -1-10420 is measured by 
the excursion from the line-center, and is 158°. This an- 
gle remained constant within the errorbars during the our 
three observing epochs, whereas the observed polarization 
increased by 1.5%. Since the excursion across the line in 
QU space became larger between 2003 and 2004, we can 
attribute this increase in continuum polarization to more 
scatterings of continuum photons by free electrons rather 
than dust particles. 

The presence of the line-effect is strong evidence that 
the electron scattering region has a geometry that deviates 
from spherical symmetry. The constant intrinsic polariza- 
tion angle indicates that the structure has remained stable. 
The nature of the Ha emitting region, and especially its ge- 
ometry, has been the topic of much debate, often centering 
on the question wh ether it is spherical o r not (see the review 
of the literature bv lDavies et al.ll2007l ). iDavies etal] (120071 ) 
use integral field spectroscopy to study refiected fight off cir- 
cumstellar dust of the emission line and show that the Ha 



line emitting region is not spherically symmetric. Here we 
confirm this result. 

It is unclear however what the polarization angle corre- 
sponds to. Assuming that the scattering is optically thin, 
the structure responsible for the scattering is oriented 
perpendicular to this at an angle of 68°. A comparison 
with the larger scale structure observed by the Hubble 
Spa ce Telescope data which tra ces the refiection nebulos- 
ity (|Humphrevs et ahll 19971 , |2002l ) is inconclusive. The angle 
does not correspond straightforwardly with either the long 
or short axis. The Hubble data reveal very clumpy struc- 
tures at sub-arcsec scales and it proves impossible to link 
the even smaller scales probed by the electron scattering, 
with the clumps or larger scale structures in the imaging. 

Finally, the increased ionization in IRC -1-10420 may be 
either explained by an increase of mass close to the star, 
due to either an infall or outfiow or due to an increase in 
temperature. 



4.2 Temperature evolution 

The photometry shows some remarkable similarities between 
the two objects. Both stars have been s ubject to studies 
of their spec tral energy dist ri bution (SED. IOudmaiier et al.l 
(1996) and iHrivnak et al] (|l989l ) for IRC -h 10420 and 
HD 179821 respectively). From fits to their SEDs, it is clear 
that the near-infrared J band photometry samples their pho- 
tospheres and not the hot dust. In both cases, this J band 
magnitude has become fainter over the past decades. At the 
same time, their V band magnitudes varied only weakly, if 
at all. 

The photometric changes in J, in particular the decrease 
in y — J, can be explained as being due to bolometric cor- 
rection effects associated with an increase in stellar tem- 
perature. Fo r IRC -1-10420, th i s is consistent with pu blished 
spectra (see IOudmaiierlll998l : iKlochkova et al.lll997l ). This 
trend is less obvious more recently. It is likely that the J 
band magnitude has reached a plateau, but it is not ap- 
parent whether this implies that the temperature increase 
has halted or not. The V—J colour is not very sensitive to 
changes in spectral type earlier than AO, so an ongoing tem- 
perature increase is not ruled out based on the photometry. 

J band data of HD 179821 is much more sparse, but 
the trend is clear. The object has become fainter in J. If 
the change in V— J colour is due to changes in effective tem- 
perature, then the star has become hotter during this pe- 
ri od, in the same way as for IRC -1-10420. For an Ay of 
2 (|Hrivnak et al.|[l989l ). and combinin g the intrinsic optical 
colou rs for supergiants as provided bv lStraizvs fc Kuriliene 
( 1981]) and intrinsic near-infrared colours due to iKoornneei 
(1983), we find that the star must have evolved to an earlier 



spectral type by about one or two spectral subtypes. 

There is some ambiguity concerning the star's temper- 
ature as derived from chemica l abundance modelling. For 
example I The venin et al.l (|2000l ) favour a lower temperature, 
consistent with the G sp ectral type. This i s lOO OK less than 
iReddv fc HrivnakI (|l999l ) and lZacs etall (|l996l ) determina- 
tions using data taken earlier. Thus, it is thus hard to draw 
any conclusions on the evolution of the stellar temperature. 
The near-infrared photometric changes may be seen as inde- 
pendent evidence that the star appears to have been evolv- 
ing towards higher temperatures. Alternatively, the temper- 
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Figure 6. A Q-U diagram of HD 179821. Each point corresponds 
to the wavelength-averaged ij-band continuum values of Q and 
U taken at different epochs. There are groupings of points in two 
separate areas of the plot, one corresponding to low PA angles 
(top) and the other to high angles (bottom). The filled circle 
shows the position of the ISP (see text). The other symbols are 
as in the previous figure. 



ature changes may be due to the formation of a pseudo- 
photosphere (Smith et al. 2004). 

The main finding in this section is that both objects 
appear to be evolving towards the blue in the HR diagram 
on human timescales. Judging from the plateauing of the 
J band photometry, IRC -1-10420 may have slowed down its 
temperature evolution, or even halted altogether. It has been 
prop o sed that the object is hitting the yellow void (de Jagerl 
1 19981 : Ide Jager etHI l200ll : iHumphrev s ct al.' 2002) and/or 
the red-edge of the bi-stability zone (jSmith ot al. 2004). In 
this case, further evolution is prevented unless a large part of 
its (pseudo-)photosphere is shredded. It will be interesting 
to follow the object to see whether a large outburst indeed 
will happen. 



4.3 Polarization Changes 

A further similarity between the two objects is the sig- 
nificant polarization variability over the past decades. 
IRC -1-10420 increased its polarization (as seen in the QU 
diagram, Fig. [3} by about 2.5%. HD 179821's polarization 
changed by more than 5%, following a straight line in QU 
space (Fig. [6|. Such changes are not uncommon in stellar 
objects, for example UX Ori variables and som e T Tauri 
stars (see e.g. lOudmaiier et al ] |200ll : lGrininlll994l ). For these 
objects, rotating clumpy circumstellar disks can result in 
large polarization changes. These occur when the clumps 
move into the line of sight of the objects, reducing the op- 
tical brightness, often considerably by many magnitudes, 
and thereby increasing the contribution of polarized light. 
However, here, in the case of both objects, the changes are 
not accompanied by large changes in the optical photom- 
etry. To make matters worse, a direct consequence of the 
UX Ori scenario, to explain the strong polarizations varia- 
tions, is that the obscuring d ust is located close to t he star 
in a compact, hot structure (|Dullemond et al.ll2003l ). This 
means that a strong near-infrared excess should be present 
and readily observable. IRC -f 10420 shows evidence for a 



significant near-infrared excess due to hot dust. In contrast, 
the near-infrared photometry of HD 179821 does not show 
any sign of hot dust close to the star and is b est explained 
by photospheric emission (jHrivnak et al.lll989l ). 

4-3.1 The electron and dust scattering of IRC +10420 

The, variable, continuum polarization of IRC -1-10420 is due 
to a combination of dust and electron scattering. The intrin- 
sic PA (10°) of the circumstellar dust polarization was de- 
rived earlier from the vector connecting the ISP and the line 
centers of the Ha emission. This PA corresponds to the short 
axis observed in the larger scale nebulosity of IRC -1-10420 
and may signify a circumstellar dusty disk structure. The in- 
trinsic angle of the electron scattering region is 158° and is 
not easily identified with the larger scale structures observed 
in the HST images. In the last observing epoch (2004), the 
line depolarization vector was much stronger than previously 
observed. Although resolution effects may have affected the 
amplitude of the excursion, the much larger continuum po- 
larization in 2004 indicates that enhanced electron scatter- 
ing is responsible for the most recent increase in polarization. 
This is supported by the fact that the continuum polariza- 
tion moved along the line vector in the QU diagram (Fig. 
[3| as well. The PA measured for the QU vector between the 
two recentmost dates is 166±3°, very close to the line polar- 
ization vectors. Such an increase is not inconsistent with the 
fact that the Ha emission EW hardly changed. Polarization 
due to (optically thin) electron scattering is more sensitive 
than (optica lly thick) Ha emission for tracing ionized mate- 
rial (see e.g. iBiorkman fc Mead3l2005l '). 

How does this compare to the longer term changes? The 
polarization evolves along a more or less straight line with 
a slope corresponding to a PA of 6±2°, which closely corre- 
sponds to the intrinsic PA of the dust scat t ering material. If 
the ISP determined value by I Jones et al.l l| 19931 ) is correct, 
then it is fair to assume that a decreasing circumstellar dust 
polarization is the main contributor to the long term vari- 
ability. A probable cause for this is a clearing dust shell. This 
can be tested by producing SED model fits to photometric 
data taken at several epochs. 

In summary, by combining the intrinsic PA of both dust 
and electron scattering and using the ISP, a relatively simple 
picture seems to emerge from the variability of the polariza- 
tion of IRC -1-10420. At first a clearing dust shell dominates 
the polarization behaviour, followed by a significant increase 
of polarization due to electron scattering. The larger elec- 
tron scattering is consistent with the observed increase in Ha 
emission, while the dust clearing can be tested with mod- 
elling the circumstellar dust. 

This global picture is probably not the final answer how- 
ever. Even the data discussed here present more questions 
already. For example, it is not clear why the point taken 
by Trammell et al. in 1991 strongly deviates from all other 
values (even one taken in the same month). The Ha de- 
polarization they observe is consistent with what we find 
ten years later, and confirms that part of the polarization is 
due to electron scattering. The large continuum value would 
indicate a violent event that led to a polarization increase 
of several per cent within weeks, however there is no other 
evidence to suggest such sudden changes in other observa- 
tions of IRC -1-10420. Similarly, it is not obvious why the 
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onset of Ha line emission in the early nineties, and its corre- 
sponding polarization of 1-1.5%, is not detected in the long 
term polarimetry, which follows a more or less straight line 
from to the seventies to the current data. A speculation is 
that a varying, rotating component of electron scattering 
was present from the onset. Given the increase in Ha emis- 
sion, it is possible that the electron scattering has become 
more optically thick. If part of the photons scatter more than 
once in their escape from the electron scattering region, a 
full or partial fli p in the position angle can result (see e.g. 
IVink et aLllioOsh . 

4.3.2 But what about HD 179821? 

The changes observed for HD 179821 are less straightfor- 
ward to interpret. The major question is how large changes 
over such an extended period can occur while, apart from 
the polarization, other observations of the star show mild 
variability at best. Firstly, we note that the ISP is hardly 
likely to vary on human timescales. This is supported by the 
results from a survey of polarization standard stars which 
reveals no variability in the obs erved polarization down to 
extremely low levels (|Clarke fc Naghizadch-Khouei 1993)- 
The Ha emission is insufficiently strong to suggest the pres- 
ence of a large electron scattering region, let alone one that 
it is variable enough to explain the observed polarization 
changes. Below we discuss the several possibilities that can 
account for the polarization and its variability observed to- 
wards HD 179821. 

4.3.2.1 Asymmetric circumstellar material Having 
been forced to rule out electron scattering as being re- 
sponsible for the variability, we now consider scattering 
by circumstellar dust. The presence of dust scattering was 
shown in the imaging polarimetry bv lKastner fc Weintraubl 
(|l995l ). The fact that our, unresolved, data show a net po- 
larization would normally indicate that part of the scat- 
tering geometry deviates from spherical symmetry. How- 
ever, this is not supported by existing imaging data of the 
circumstel lar material of the object. The scattering dust 
traced by iKastner fc Weintraubl (119951 ) stretches over sev- 
eral arcseconds, but appears to be distr ibuted spherically 
symmetrically. HST images published bv lUeta et al.l (|2000l ) 
show small scale structure but largely a round appearance. 
iBuiarrabal et all 's Il992l CO rotational data maps show a 
slightly elongated structure at larger (~5 arcsec) distances 
from the star at a PA of 50° . 

Even if the slight asymmetry would be responsible for 
the observed net polarization, it is hard to reconcile this 
with the observed variabi lity. With a kinematic age of order 
thousands of years (e.g. iKastner fc Weintraubl [l995l '). one 
would not expect changes at a fraction of this period, i.e. 
15 years, to still be traceable. A more practical objection to 
the extended emission being variable is that the intensity of 
the scattered light rapidly drops as a function of distance 
to the star, so there is simply not enough light available to 
induce variations of 4% in polarization in the total light. 

A final possibility is that a compact asymmetric struc- 
ture much closer to the star is responsible for the variable po- 
larization. As mentioned earlier, the SED of the object leaves 
no room for a significant, hot dusty component. A caveat is 



of course that the near-infrared photometry of the object is 
limited, so any sudden, short-lived events could have been 
easily missed. If the ISP value is, as the data suggest, in 
the upper half of the QU diagram, the major polarization 
activity occurs when the object occupies the lower half of 
the QU diagram. A prediction would then be that during 
such periods, excess radiation in the near-infrared due to 
hot dust would be present. NIR data of HD 179821 is very 
sparse, but the last set of J, K photometry (from 2MASS), 
was taken during Melikian's observing period. These data 
follow the gradual faintening of the star in the NIR and do 
not indicate any excess emission. In addition, the optical 

photometry appears steady throughout. 

F or completeness, we do note that lOudmaiier et al] 
a995l ) found variable first-overtone CO emission at 2.3 /xm 
towards this object. This emission typically arises in hot, 
dense material within the dust sublimation radius. Oudmai- 
jer et al. interpreted this as variable mass loss. If the CO 
emission originates from a geometry that is not spherically 
symmetric, it could be associated with ejections of dense 
gas similar to, but shorter than, for example the asymmet- 
ric m ass loss event observed towards HD 45677 ( Patel et al] 
l2006h . The resulting situation could give rise to Rayleigh 
scattering and a net polarization of light. However, whether 
this would give rise to the large polarization variations is not 
clear. We do point out that CO emission from the yellow hy- 
pergiant p Cas could also be expla ined by shocks due t o pul- 
sations in the stellar photosphere (ICorlova et al I2OO6I ). This 
brings us to the central star as the origin of the variations. 

4.3.2.2 An asymmetric star illuminating the dust 
shell Above, we considered the case of a central source 
surrounded by asymmetric scattering media. The asymme- 
tries involved result in a net polarization of the total light. 
However, we tacitly assumed the central object to be either 
a point source or an isotropically emitting sphere. 

It is well known that when a central source illuminates 
a spherically symmetric envelope, any asymmetries in the 
central sour ce can induc e a net polarization. As for example 
reviewed bv lBovle et al.l 1 )1986 ), who find that a photosphere 
with a nonuniform surface brightness can result in net po- 
larization when its radiation is scattered off circumstellar 
dust. Other examples are changes found in the polarization 
spectra of evolved stars where it has been inferred that dif- 
ferent molec ular bands cover di fferent parts of the stellar 
photosphere (|Bieging et al.|[2"006l ). The polarization changes 
over molecular bands can be explained by both Rayleigh 
scattering by particles in clouds high in the atmosphere and 
scattering due to circumstellar dust (|Raveendranlll99ll ). In 
all these cases the break from symmetry is not due to the cir- 
cumstellar material, but due to the asymmetric illuminating 
source. 

The extensive photometric observations and period 
analysis of lLe CoroUer et al.l (120031 ) indicate that HD 179821 
is a non-radial pulsator (NRP). It is not yet clear whether 
the shape of the star changes altogether, as we lack the spec- 
tral data (Le CoroUer et al. quote unpublished data implying 
that the velocity variations are modest), or due to changes in 
surface temperature, the situation Le CoroUer et al. prefer. 
In either case, the circumstellar material will be irradiated 
by an asymmetric source. Then, the polarization changes can 
be expected to occur around a preferred axis, as, apart from 
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rotation, the pulsation nodes are not expected to stochasti- 
cally change position on the star. 

In order to investigate whether it is plausible that a 
star with an aspherical shape can yield observable polariza- 
tion, we performed some test calculations of a star covered 
by large spots, surrounded by a spherically symmetric dust 
shell. The model set-up consists of a star with a spot. In the 
situation where one (bright) spot is present, the polarization 
will increase with viewing angle. That is to say, if the spot 
is located in the line of sight, we will see its light directly 
and not much net polarization will be measured. If on the 
other hand the system is inclined, less direct light will be 
seen and the polarized light contribution increases, hence 
a larger polarization is observed. We should also note that 
unless the scattering occurs very close to the star (at a few 
stellar radii), the polarization percentage is independent of 
the radius of the shell. This is due to the constant geometry 
of the situation. 

For the simple case where all the stellar flux is radi- 
ated from one spot covering 10 per cent of the stellar sur- 
face, we flnd that polarizations up to 3 per cent are readily 
produced when the star is mildly inclined. For smaller sur- 
face areas, the polarization can be even larger, exceeding 
10 per cent. This result c onfirms earlier studies which used 
less e xtreme parameters (|Raveendranlll99l] : lAl-Malki et al.l 
1 19991 ). Of course this is an idealized situation, but it demon- 
strates that net polarizations of order per cent are certainly 
possible. Indeed, from an observational point of view, po- 
larization changes of around 1 pe r cent have been observed 
towards Betelgeuse (|Haveslll98A . These changes were due 
to variable hotspots on the stellar surface. Therefore, polar- 
ization changes can indeed be of order few per cent if the 
star itself is anisotropic. 

One could argue why the polarization variability is not 
modulated by the 140 or 200 day period. Perhaps it is on 
the shorter term, but we lack the data to properly inves- 
tigate this. However, where we could match photometric 
and polarimetric data, on shorter timescales it appears that 
the polarization values closely follow the photometry. It is 
therefore possible that the polarization is variable at simi- 
lar timescales to the pulsations. A further question is why 
such long term changes are observed. We note that the mean 
photometry shows long term changes, probably tracing the 
global evolution of the photosphere, which would have an ef- 
fect on the non-radial pulsations. Significantly, the preferred 
plane would not necessarily change, and this is found. 

A critical test of the scenarios discussed above is si- 
multaneous photo-polarimetric monitoring of HD 179821. If 
non-radial pulsations are responsible for the observed po- 
larization, both observables should be strongly correlated. 
If, on the other hand, episodic mass ejections are the main 
cause of the variations, we would expect enhanced near- 
infrared emission due to hot dust during significant polar- 
ization changes. 



5 CONCLUSION 

We have presented spectropolarimetry and long term photo- 
polarimetry for two post-Red Supergiants. A strong depolar- 
ization across the Ha emission line is found for IRC -1-10420, 
suggesting an electron-scattering region that is not circularly 



symmetric, confirming the results of iDavies et al.l |20o3), 
who found such evidence from their integral field spec- 
troscopy of the object. The time evolution indicates that 
the source has increased in temperature until at least 1995, 
after which the J band photometry indicates this may have 
levelled off. If the temperature increase of the object has in- 
deed halted, an increase in mass loss or mass infall rate can 
be responsible for the observed Hq emission and polariza- 
tion. 

HD 179821, a cooler object, has less Ha emission, and 
no depolarization across the Ha line could be detected. The 
photometry implies that this evolved object is also under- 
going a change in temperature. Hitherto, temperature de- 
terminations of the object were inconclusive, but if the star 
is a G supergiant, then the observed change in photome- 
try suggests it has become earlier by one or two subclasses. 
Strong changes at the 5 per cent level in polarization over 
the past 15 years have been detected. During the same time, 
the optical photometry has only varied by at most 0.2 mag- 
nitudes. The most obvious explanations for this observation 
such as changes in either electron or dust scattering can be 
ruled out. A complication is that the polarization is not cor- 
related with any other observable. The possibility that the 
star undergoes irregular, asymmetric, mass ejections is dis- 
cussed. However, there is little evidence for the occurrence 
of such events. Instead, it is proposed that the star itself 
is asymmetric and that its anisotropic radiation, which is 
scattered off the circumstellar dust results in the net polar- 
ization observed. 
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